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Abstract. We find chimera states with respect to amplitude dynamics in a network of Stuart- 
Landau oscillators. These partially coherent and partially incoherent spatio-temporal patterns 
appear due to the interplay of nonlocal network topology and symmetry-breaking coupling. As 
the coupling range is increased, the oscillations are quenched, amplitude chimeras disappear 
and the network enters a symmetry-breaking stationary state. This particular regime is a 
novel pattern which we call chimera death. It is characterized by the coexistence of spatially 
coherent and incoherent inhomogeneous steady states and therefore combines the features of 
chimera state and oscillation death. Additionally, we show two different transition scenarios 
from amplitude chimera to chimera death. Moreover, for amplitude chimeras we uncover the 
mechanism of transition towards in-phase synchronized regime and discuss the role of initial 
conditions. 


1. Introduction 

Chimera states represent an intriguing phenomenon in a dynamical network of coupled elements, 
which spontaneously separates into two coexisting domains with dramatically different behavior 
[1]. They occur, surprisingly, in networks of identical units and symmetric coupling schemes 
mm- In spite of the fact that these symmetry-breaking states have been the subject of many 
investigations, e.g., msiEiiziiHiiaiiniiiiiiiziiiainiiis], there still remain many issues to 
be addressed. A recent review on chimera states m discusses a number of significant open 
questions concerning, for instance, the necessary conditions for their existence or generalizations 
of chimera states. A generalization to amplitude dynamics and to steady state coherence- 
incoherence patterns has been provided in [1]. In particular, in [I] it has been shown that 
chimeras of steady states appear due to a special symmetry-breaking coupling which originates 
from oscillation death investigations [m dsi ng EO]. Oscillation death, i.e., the stabilization 
of inhomogeneous steady states in coupled oscillator networks, appears when a homogeneous 
steady state splits into at least two distinct symmetry-breaking branches - upper and lower m- 
In comparison with amplitude death, which is another type of oscillation suppression, i.e., the 
stabilization of an unstable homogeneous steady state, oscillation death is much less investigated 

m 

Initially found for phase oscillators [2], chimera states can be treated as long chaotic transients 
towards the in-phase synchronized regime. In this context the network size plays an essential 
role. In particular, it has been shown theoretically m and confirmed experimentally |23| that 
the chimera lifetime grows exponentially with the system size. Therefore, chimera states are 


stable for infinitely large networks and transient for finite-size ensembles. Recently, in |23] the 
question of what is the smallest possible number of oscillators allowing for chimera states has 
been addressed. Another important problem, also considered as a necessary condition for the 
existence of chimera states, is the choice of initial conditions. Random initial conditions do not 
always guarantee chimera behavior. This is due to the fact that classical chimera states typically 
coexist with the completely synchronized regime. The coexistence of two stable solutions is a 
signature of bistability which is typically associated with hysteresis. The basin of attraction for 
chimera states can be relatively small in comparison with that of the synchronized state. For 
that reason chimera states remained for a long time undetected. Indeed, the first experimental 
studies pHl 071 IHni ITT] on chimera states were provided only twelve years after their 

theoretical discovery. Consequently, it is reasonable to use specially prepared initial conditions 
to ensure chimera patterns. However, it is worth noting that a chimera state is not just a 
temporary trace of initial conditions which disappears in time, but a persisting pattern with a 
long lifetime. 

In the present work we investigate the connection between two symmetry-breaking effects: 
chimera states and oscillation death. We show that the interplay of nonlocal coupling topology 
and symmetry-breaking coupling leads to a number of novel partially coherent inhomogeneous 
spatial patterns, for instance amplitude chimeras and chimera death. Additionally, we analyze 
the transition scenarios between these two patterns. Moreover, the transient amplitude chimeras 
and the transition mechanism from amplitude chimera to a synchronized regime is investigated. 
Further, the newly found amplitude chimeras are analyzed with respect to initial conditions and 
transients. In particular, we aim to understand the mechanism of transition from amplitude 
chimera to in-phase synchronization. Finally, we address the question of how the lifetime of 
amplitude chimeras depends on initial conditions and what are the optimal initial distributions. 

2. Model 

We analyze the paradigmatic model of Stuart-Landau oscillators [21 [3211331 [311 ESI [361137]: 

i = f{z) = {X-\- iuj — (1) 

where z = re'^^ = x -h i?/ G C, A, cj G R. For A > 0, the uncoupled system exhibits self-sustained 
limit cycle oscillations with radius tq = and frequency oj. Therefore, the Stuart-Landau 
system represents a generic model for nonlinear oscillators close to a Hopf bifurcation. We 
investigate a ring of N nonlocally coupled Stuart-Landau oscillators [1]: 

(j 

= + ^ (Rezk-Rezj), ( 2 ) 

k=j-P 

where j = 1, 2,..., A^. The coupling parameters, which are identical for all links, are the coupling 
strength cr G R and the coupling range P/A, where P corresponds to the number of nearest 
neighbors in each direction on a ring. Here we consider coupling only in the real parts, since 
this breaks the rotational symmetry of the system which is a necessary condition for the 
existence of nontrivial steady states Zj ^ 0 and thus for oscillation death m- With respect 
to applications this means that the oscillators are coupled only through a single real variable 
X. If compared with slow-fast systems such as, for example, the FitzHugh-Nagumo oscillator, 
which demonstrates two separated time scales, the Stuart-Landau oscillator is a generic model 
for sinusoidal oscillations which have only one characteristic time scale defined by the frequency 


3. Amplitude chimeras and chimera death states 

While tuning the coupling range P for a fixed value of the coupling strength a we uncover a 
variety of dynamical regimes in equation (2) which are shown as space-time plots, color-coded 
by the variable ?/, and snapshots in Fig. In particular, we find chimera behavior with respect 
to amplitude dynamics, i.e., amplitude chimeras [1], where one part of the network is oscillating 
with spatially coherent amplitude, while the other displays oscillations with spatially incoherent 
amplitudes (Fig. & ,f). It is important to note that within the incoherent domain of the 
amplitude chimera the center of mass of each oscillator Zc.m. = Jq ^jit)dt/T^ where T — 27r/uj 
is the oscillation period, is shifted away from the origin, while the elements from the coherent 
subgroup oscillate around the origin. To illustrate this fact we calculate the distance between 
the center of mass of each oscillator and the origin rc.m. (Fig. [^a). Surprisingly, the shape 
of this profile is very similar to that of the mean phase velocity profile for the classical phase 
chimeras. However, for amplitude chimeras the mean phase velocity profile Oj is flat (Fig. If)- 
This fact discloses the crucial feature of amplitude chimera: averaged phase velocities remain 
the same for every element of the network, since the phases are correlated even within the 
incoherent domain. In more detail. Fig. shows the phase portraits of all oscillators in the 
complex 2 ; plane clearly demonstrating the limit cycles with different amplitudes and centers of 
mass. The nodes from the coherent part oscillating with larger amplitudes around the origin 
perform fast motions in the complex z = x-\-iy phase plane, while the elements of the incoherent 
domain with smaller amplitudes are slowed down. Therefore, the angular frequencies of all the 
nodes on average remain the same. Consequently, we observe pure amplitude chimera - chimera 
behavior exclusively with respect to amplitude dynamics rather than the phase, in contrast 
to amplitude-mediated chimeras, for which both phase and amplitude are in a chimera state 
[291 EE]. 

The increase of the coupling range P/N induces a transition from amplitude chimeras to 
an in-phase synchronized state (Fig. ,g). By increasing the coupling range even further we 
detect a novel pattern which provides bridging between chimera states and oscillation death. 
Therefore, we call it chimera death P (Fig- -j). In this regime the oscillations die out in 

a peculiar way. The population of identical oscillators breaks up into two domains: (i) spatially 
coherent oscillation death, where the neighboring elements of the network are correlated forming 
a regular inhomogeneous steady state, and (ii) spatially incoherent oscillation death, where the 
sequence of populated branches of the inhomogeneous steady state of neighboring nodes is 
completely random. It is important to note that the term “coherent/incoherent” refers to the 
coherence/incoherence in space, i.e., spatial correlation, which should be distinguished from 
temporal coherence which refers to correlations in time of the dynamics. Because of symmetry 
reasons, for a node on the upper branch -hi of the inhomogeneous steady state in the 

left half of the system, there always exists a mirror state shifted by phase tt (anti-phase) in 
the complex plane, i.e., located on the lower branch ~ —1, in the right half of the system. 
Interestingly, the increase of the coupling range for fixed coupling strength also induces structural 
changes of the chimera death pattern: the coherent spatial domain may consist of one or several 
clusters where the nodes are on the same branch of the inhomogeneous steady state. With 
increasing coupling range the number of clusters in the coherent spatial domain is decreased, 
see the scenario in Fig. from (c,h) (5 clusters) via (d,i) (3 clusters) to (e,j) (1 cluster). The 
illustration of the chimera death pattern with the maximum number of clusters in the coherent 
part is provided by Fig. [^-c. The phase portrait typical for chimera death is shown in Fig. [^, 
it consists of two fixed points shifted by a phase tt, corresponding to the two branches 
of the inhomogeneous steady state. It should be noted that these two fixed points are very close 
to the centers of mass of those oscillations in Fig. with the smallest amplitude, i.e., in the 
center of the incoherent domain of the amplitude chimera. 

To provide an overall view on the network behavior for the wide range of coupling parameters 


we plot the map of regimes in the plane of coupling range and coupling strength. The oscillatory 
behavior of the network is represented by amplitude chimeras (blue region in Fig. Q, which 
are observed for small coupling range, and by in-phase synchronized oscillations (light green 
region in Fig. |^. The steady state solutions occur for larger values of coupling parameters and 
are manifested by chimera death (red regions in Fig. with different hatching). The region 
of chimera death is divided into several regimes depending on the number of coherent clusters. 
Chimera death may consist of two coherent domains, which are formed by one cluster on the 
upper and one on the lower branch of the steady state (dark and light red stripes in Fig. [^. 
Moreover, the coherent regions may split into three clusters each (dark red and yellow stripes 
in Fig. ??) or more (dark red in Fig. [^. The typical space-time patterns and snapshots in 
panels (c,d,e) and (h,i,j) of Fig [^provide illustrations of the different multi-cluster chimera death 
states, and the corresponding coupling parameters are indicated in the phase diagram (Fig[^ 
by empty diamonds. It is important to note that the chimera death regime is characterized by 
high multistability. Therefore, chimera death patterns with different number of clusters in the 
coherent part coexist. The map of regimes strongly depends on initial conditions, in particular, 
the borders between chimera death patterns with different number of clusters. However, the 
boundary separating the oscillatory regime from the steady state regimes (with multiple cluster 
numbers) appears to be less sensitive to initial conditions and remains almost the same for 
different realizations of initial conditions. The existence of two distinct transition scenarios 
from the oscillatory to the steady state regime becomes evident from Fig. For a small 
value of coupling strength cr = 10 the amplitude chimera gives way to the chimera death state 
after passing through in-phase synchronized oscillations when the coupling range is increased 
(diamonds in Fig. Q. In contrast, for a large value of the coupling strength, for example, 
a = 26 a slight increase of the coupling range from P/N = 0.04 to P/N = 0.05 destroys 
amplitude chimeras and directly leads to chimera death. 

4. Transient behavior of amplitude chimeras 

An important feature of classical chimera states is their transient character and strong 
dependence on initial conditions. Our numerical analysis shows that amplitude chimeras 
although manifesting a long lifetime (up to 10^ dimensionless time units) for the optimal choice 
of initial conditions, can rapidly transform into in-phase synchronized regime for completely 
random initial conditions. 

To gain more insight into the transient behavior of amplitude chimeras we trace the system’s 
evolution in time having fixed all other parameters. Figure [^-d shows corresponding snapshots 
for different time values. As initial condition here we take chimera death (t = 0). The pattern 
at an early stage {t = 4) shown in Fig. is reminiscent of the initial conditions, while the 
regime of amplitude chimeras (Fig. [^) is reached later {t = 65). However, after approximately 
125 time steps it starts deforming (Fig. until it turns into an in-phase synchronized pattern 
(Fig. [5ji). This is also visible in the space-time plot (Fig. [^). 

To disclose the transition mechanism from amplitude chimera to the in-phase synchronized 
regime we analyze phase portraits in the amplitude chimera state (Fig. [^) and during the 
transition (Fig. ^jp) for three selected nodes of the network: one from the coherent part (j = 50) 
and the other two from the incoherent domain {j = 99,102). The time ranges covered in Fig. 
and b, respectively, are indicated by two green vertical lines in the space-time plot of Fig. [^. In 
the regime of amplitude chimera the node from the coherent part (j = 50) oscillates with large 
amplitude around the origin, while the two nodes in the incoherent domain (j = 99,102) oscillate 
with much smaller amplitudes around their own centers of mass shifted from the origin in two 
opposite directions (Fig. |§i). During the transition the amplitudes of the incoherent oscillations 
grow till they reach the value of the coherent oscillations. At the same time the centers of mass 
for nodes from the incoherent domain move towards the origin (Fig. |§)). Therefore, during the 



Figure 1. Space-time plots (left panel) and snapshots for the variable yj{t) (right panel) 
in a network of Stuart-Landau oscillators for coupling strength a = 10 and varying nearest 
neighbors number P. (a),(f) P = 4: amplitude chimera; (b),(g) P — 10: in-phase synchronized 
oscillations; (c),(h) P = 15: multi-cluster (> 3) chimera death; (d),(i) P = 25: 3-cluster chimera 
death; (e),(j) P = 45: 1-cluster chimera death. Other parameters: N = 100, A = 1, cj = 2. 



Figure 2. Amplitude chimera: (a) snapshot at t = 1000 for rc.m. (the distance between the 
center of mass for every oscillator and the origin); (b) mean phase velocity profile; (e) phase 
portraits of all oscillators in the complex z — x-\-iy plane. Parameters: N — 100, P = 4, a = 14, 
A = 1, cj = 2. 
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Figure 3. Multi-cluster (> 3) chimera death: (a) space-time plot for the variable yj] (b) 
snapshot for the variable yj] (c) phase portrait of all oscillators in the complex z — x-\-iy plane. 
Parameters: N — 100, P = 5, a = 26, A = 1, cj = 2. 



Figure 4. Map of dynamical regimes for N = 100, A = 1, cj = 2 in the plane of coupling 
range P/N and coupling strength a for specially prepared initial conditions, showing 1-cluster 
chimera death (1-CD); 3-cluster chimera death (3-CD); multi-cluster (> 3) chimera death {N- 
CD); amplitude chimera (AC); in-phase synchronized oscillations (SYNC). Diamonds mark the 
parameter values chosen in Fig. Other parameters: N = 100, A = 1, cj = 2. The detailed 
explanation of initial conditions is given in Section 5. 


transition the phase trajectories for the nodes j — 99,102 have spiral shapes (Fig. |§)) . 

An important observation is that the transition mechanism described above depends on the 
initial conditions. For example, asymmetric initial conditions may lead to another scenario 
which involves an asymmetric state. In this case nodes from the incoherent part join one by one 
their neighboring synchronized state. 

5. The role of initial conditions 

The first initial condition we choose originates from the oscillation death problem and is typically 
used for the inhomogeneous steady state m- Specifically, the system is divided into two equal 
domains: one is located on the upper branch and the other is on the lower branch. In our 





















































particular example of = 200 nodes this implies that one half of the network (j = 1,100) is 
on the upper branch ^ +1 while the other half is on the lower branch ^ —1 (Fig. [^). 

In this case the lifetime of amplitude chimeras is t < 200 and decreases significantly when some 
small random shift from the branches is applied for every node. 

Another type of initial condition we have tested is the modification of the previously discussed 
set by adding two incoherent regions, where the nodes are distributed randomly (Fig. lit) . This 
leads to dramatically short lifetime (t < 30) of the amplitude chimera, which rapidly transforms 
into an in-phase synchronized state or a traveling wave. 

Next we allow randomness in initial conditions for all nodes of the network, however, observing 
certain symmetries. In particular, in Fig. the positions for 1/4 of the nodes j = 1,...,50 
are chosen randomly around the upper branch, then these values are mirrored to the nodes 
j = 51,...,100 by setting zj — Finally, the positions of the other half of the nodes 

j = 101, ...,200 are obtained by applying a phase shift of tt which we call “anti-phase partner” 
condition: zj = —^j+Ar/2 with j mod N. The initial condition scheme shown in Fig. [71 is 
constructed in a similar way as in Fig. [^. The only difference is that the values for thenrst 
j = 1,...,50 nodes are chosen randomly around both branches of the inhomogeneous steady 
state (and not only around the upper branch like in Fig. 14)- This amendment does not 
have any impact on the lifetime of amplitude chimeras, which is on the other hand strongly 
affected by the p resence of symmetries in the initial conditions. Our simulations suggest that 
amplitude chimeras tend to have very short lifetime if in the incoherent domain the initial 
positions of nodes near the upper branch are uncorrelated with those close to the lower branch. 
An increase of the amplitude chimera lifetime can be achieved by introducing symmetries in the 
initial distribution (“anti-phase partner” condition). Therefore, optimal results with t > 3000 
are obtained for random initial conditions with symmetries (Fig. [^,d). In conclusion, our 
numerical analysis shows that for random distributions without symmetries amplitude chimeras 
appear to be short transients towards in-phase synchronized regime, while symmetric conditions 
significantly increase their lifetime up to t > 10^. Moreover, random initial distributions may 
also lead to asymmetric spatio-temporal patterns. 

To provide a comprehensive view on the impact of initial conditions we compare the map of 
regimes calculated for a specially prepared initial set (Fig. and for random initial conditions 
(Fig. [^. For the diagram shown in Fig. we use as initial condition an amplitude chimera 
profile at a fixed time (similar to Fig. [^), which is obtained in the following way: For a fixed set 
of parameters in the amplitude chimera regime, the system is divided into two equal domains; 
half of the nodes are located on the upper branch and the other half is on the lower branch 
of the inhomogeneous steady state (similar to Fig. [^); this initial condition then evolves into 
an amplitude chimera state which is used as initial condition for all other parameter values 
{(j^P/N). The random initial distribution (similar to Fig. [^) leads to a significantly different 
phase diagram (Fig. in comparison with the specially prepared initial conditions (Fig. Q. 
The domain of the amplitude chimera is decreased for strong coupling if compared with Fig. 

Moreover, the regime of chimera death obtained from carefully chosen initial conditions 
(Fig. H) is replaced by the regime of incoherent oscillation death (inset in Fig. which results 
from random initial conditions. It should be noted that the map of regimes shown in Fig. 
is obtained for one single realization of random initial distribution and no averaging of initial 
conditions has been provided. Therefore, it is not statistically representative that chimera death 
is not observed in Fig. [^ This is merely due to the fact that the probability to achieve chimera 
patterns starting from one random realization of the initial conditio n is very low. However, the 
occurrence of chimeras for random initial conditions, in general, is not excluded. 




Figure 5. Transient amplitude chimera for N = 200, A = l,6(; = 2, P = 5 and a = 20. 
Snapshots for the variable yj (left panel): (a) t = 160: in-phase synchronized oscillations; (b) 
t = 140: deformed amplitude chimera at the transition point; (c) t = 65: amplitude chimera; 
(d) t = 4: trace of initial condition. The right panel shows the space-time plot for the variable 
yj{t). Blue diamonds mark the time values chosen for the snapshots in the left panel. The green 
lines mark the time ranges chosen for the phase portraits in Fig. 


6. Conclusions 

In conclusion, we have provided a connection between two symmetry-breaking effects - chimera 
states and oscillation death. In particular, we uncover new spatio-temporal patterns, i.e., 
amplitude chimeras and chimera death. Amplitude chimeras represent a generalization of 
chimera behavior to amplitude dynamics, and chimera death generalizes chimeras to steady 
states. It is shown that different transition scenarios from amplitude chimeras to chimera death 
are possible. Chimera death patterns with different numbers of clusters in the coherent part 
exhibit a high degree of multistability, which is related to a variety of hysteretic scenarios. 
Moreover, the transient dynamics of amplitude chimera is investigated disclosing the transition 
mechanism from amplitude chimera to the in-phase synchronized regime. Additionally, we 
find optimal initial conditions for amplitude chimeras showing that random distributions with 
particular symmetries essentially enlarge their lifetime. 
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Figure 6. Phase portraits of three selected nodes of the network in the complex z = x + iy plane 
for (a) amplitude chimera state and (b) during the transition to in-phase synchronization. Red 
color marks the trajectory of node j = 50 (from the coherent part). Dark and light blue colors 
show trajectories of nodes from the incoherent domain j = 99 and j = 102 correspondingly. 
Time ranges as marked in green in Fig. Parameters: N = 200, A = 1, cj = 2, P = 5, cr = 20. 
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Figure 7. Initial conditions for amplitude chimeras used for N = 200, A = l,6(; = 2, P = 4 and 
a — 20: (a) coherent distribution: one half of the nodes on the upper and one half on the lower 
branch of the inhomogeneous steady state (amplitude chimera lifetime t — 130); (b) combination 
of coherent and incoherent domains with no symmetries in the randomly defined incoherent part 
(amplitude chimera lifetime t = 25); (c),(d) random initial conditions with symmetries (optimal 
amplitude chimera lifetime t > 3000). 
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Figure 8. Map of dynamical regimes for N = 100, A = 1, cj = 2 in the plane of coupling range 
P/N and coupling strength a for random initial conditions as in Fig. [^d): amplitude chimera 
(AC); in-phase synchronized oscillations (SYNC); incoherent oscillation death (Incoherent OD). 
The inset shows a space-time plot for the variable yj{t) for the coupling strength cr = 30 and 
coupling range P/N = 0.5. 
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